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Abstract 

Ordering ambiguity associated with the von Roos position dependent 
mass (PDM) Hamiltonian is considered. An affine locally scaled first order 
differential introduced, in Eq.(9), as a PDM-pseudo- momentum operator. 
Upon intertwining our Hamiltonian, which is the sum of the square of this 
operator and the potential function, with the von Roos d-dimensional 
PDM-Hamiltonian, we observed that the so-called von Roos ambiguity 
parameters are strictly determined, but not necessarily unique. Our new 
ambiguity parameters' setting is subjected to Dutra's and Almeida's [11] 
reliability test and classified as good ordering. 

PACS numbers: 03.65.Ge, 03.65.Fd,03.65.Ca 

1 Introduction 

Quantum mechanical Hamiltonians with position dependent mass (PDM) con- 
stitute interesting and useful models for the study of many physical problems 
[1-31]. They are used in the energy density many-body problem [1], in the de- 
termination of the electronic properties of the semiconductors [2] and quantum 
dots [3], in quantum liquids [4], in 3 He clusters [5] and metal clusters [6], in 
the Bohmian approach to quantum theory (cf., e.g. [7]), in the full and partial 
wave-packet revivals(cf., e.g., [8]), etc. Comprehensive reviews on the applica- 
bility of such position dependent mass settings could be found in the sample of 
references in [1-12] 

However, it is concreted that an ordering ambiguity conflict arises in the 
process of defining a unique kinetic energy operator, due non-commutativity 
between the momentum operator p x — —id x and the position dependent mass 
M (x) — m m(x). A problem that has shown poor advancement over the last 
few decades. 
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In general, working on PDM Hamiltonians is inspired by the von Roos Hamil- 
tonian [15] proposal (with ft = 2m = 1) 
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m (x) a d x m (x) 13 d x m (x) 7 + m (x) 7 d x m (x) 13 d x m (x) a +V(x) (1) 

subjected to von Roos constraint 

a + p + ~/= -1 ; tt,ft 7 el (2) 

Hamiltonian (1) may, in a straightforward manner and with the constraint (2), 
be very well recast ( cf., e.g., [10]) as 

H=-d x (^—]d x + V(x), (3) 



m (x) 



with 



f 

2 m (x) m (x) 



V (x) = ^(1 + [a(a + (3 + 1) + (3 + 1]^J + V (x) , (4) 



where primes denote derivatives. Obviously, nevertheless, the profile of V (x) 
(namely the first two terms in (4)) changes as the parameters a, (3, and 7 
change, manifesting therefore the eruption of ordering ambiguity in the process 
of choosing the kinetic energy operator T. Hence, a, (3, and 7 are usually called 
the von Roos ambiguity parameters. 

Several proposals for the kinetic energy operator are suggested in literature. 
Amongst exist; the Gora and Williams (/3 = 7 = 0, a = — 1) [16], Ben Danial 
and Duke (a = 7 = 0, (3 = -1) [17], Zhu and Krocmcr (a = 7 = -1/2, (3 = 0) 
[18], and Li and Kuhn (/? = 7 = —1/2, a = 0) [19]. However, the Hermiticity of 
the kinetic energy operator, the current density conservation, the experimental 
results [20-25], and the condensed matter theories [26,27] may give some ideas on 
the identity of the von Roos ambiguity parameters. Applying Hamiltonian (1) 
to an abrupt heterojunction between two crystals (cf., e.g., sample of references 
in [25-28,30]), for example, implied that for a ^ 7 the wave function vanishes 
at the heterojunction (i.e., the heterojunction plays the role of an impenetrable 
barrier). Hence, the only feasible cases are due a = 7 to ensure the continuity 
of m (x) a %j) (x) and m (x) a+ [d x ip (x)] at the heterojunction boundary. 

Very recently, however, Dutra and Almeida [11] have carried out a reliability 
test on the orderings available in literature. They have used an exactly solvable 
Morse model and concluded that the orderings of Gora and Williams (a = (3 = 
7 = 0, a = -1) [16], and Ben Danial and Duke (a = a = 7 = 0, = -1) [17] 
should be discarded for they result in complex energies. Nevertheless, they have 
classified the ordering of Zhu and Kroemcr (a — 0, a = 7 = —1/2, (3 = 0) [18], 
and that of Li and Kuhn (a ^ a ~ 0, /3 = 7 = —1/2) [19] as good orderings. 
Yet, they have shown that Weyl (cf., e.g., Borges in [30]) and Li and Kuhn [19] 
orderings are equivalent. 
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Ultimately, therefore, the continuity conditions at the heterojunction bound- 
aries and Dutra's and Almeida's [11] reliability test single out Zhu and Kroemer 
(a = 0, a = 7 = — 1/2, /3 = 0) [18] as good ordering. This, in effect, inspires our 
current methodical proposal and manifests the introduction of a PDM-pseudo- 
momentum operator which, in turn, leads to a new good ordering. 

On the other hand, within a Liouvillean-type change of variables spiritual 
lines, the point canonical transformation (PCT) method for Schrodinger equa- 
tion often mediates a transition between two different effective potentials. That 
is, in the PCT settings, one needs the exact solution of a potential model in a 
class of shape invariant potentials to form the so-called reference/old potential. 
The reference/old potential along with its exact solution (i.e. eigenvalues and 
eigenfunctions) is then mapped into the so-called target/new potential , hence 
exact solution for the target/new potential is obtained. For more details on 
this issue the reader may refer to , e.g., ref. [14]. 

In this parer, we recollect (in section 2) the d-dimensional form of the von 
Roos Hamiltonian suggested by Quesne [10] and introduce our PDM-pseudo- 
momentum operator. The consequences of such operator's setting on the von 
Roos ambiguity parameters are given in the same section. Moreover, the corre- 
sponding d-dimensional radial Schrodinger Hamiltonian and the PCT d-dimcnsional 
mapping are also reported. Our concluding remarks are given in section 3. 

2 (/-dimensional von Roos Hamiltonian and PDM- 
pseudo-momentum operators 

Quesne in [10] has suggested a general form of von Roos e?-dimensional PDM 
Schrodinger equation 

'"m (q) Q dj-m (q)' 3 djin (q) 7 + m (q) 7 djtn (q) 13 djin (q) Q | V (q) 

+ 0/(q)-£}V(q) = 0, (5) 



where q = ((71,(72,' •• ,1d), dj = d/dqj, j = 1,2, ••• ,d, m (q) is the dimen- 
sionless form of the mass M (q) = m m (q), V (q) is the potential function, and 
summation runs over repeated indices. In this case, the (f-dimcnsional PDM 
Schrodinger Hamiltonian reads 

with 

V(q) = i(l + / 3 )^^-[a(a + / 3+ l) + / 3+ l]t^^ + V( q) .(7) 
z m (q) m (q) 

Let us now consider, for simplicity, quasi-free-particles' setting (i.e., V (q) = 
0). Then it would be obvious that the quasi- free-particles' Hamiltonian structure 
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suggests that the kinetic energy operator 

\m(qjy 2 m(q) m (q) 

(8) 

may, mathematically speaking, very well be expressed as the square of a first- 
order differential vector operator of a general form 

il j = -i{F(m(d))d j +G j (m(ci))}. (9) 

This would (with F(m(q)) = F (q) = F, (m(q)) = Gj (q) = Gj for sim- 
plicity, and T = ft 2 = ^jjiTjlTj) imply 

f = —F 2 djdj— [F (djF) + 2FGj] d 3 - [F {djGj) + GjGj] , (10) 
If we compare Eq.(10) with (8) we obtain 

F(q) = ±- 



V m (q) ' 

= _2F djF = _ [F {djF) + 2FGj] , 

m(q) 

= ^ ; (11) 

The structure of our first-order differentioal operator is therefore clear and 
can be cast as 

tl j = -i^F(d)d j + ^[d j F(c l )}Y (12) 

At this point, it should be noted that our operator II is Hermitian and represents 
the position-dependent-mass generalization of the ordinary momentum operator 
Pj = —idj (i.e., at constant mass settings M (x) — m ). Hence, II could be 
labeled, hereinafter, as a PDM-pseudo-momentum operator. 

2.1 Consequences of our PDM-pseudo-momentum oper- 
ator II on the von Roos ambiguity parameters 

In a straightforward manner it is easy to show that 

1 djdjmjq) _ 7 [djm(q) ] 2 
m{q)J " J ' 4 m(q) 2 16 m(q) 3 



ft 2 = - d . * 0. + - T7, 3 1,™ ■ ( 13 ) 



Comparing this result with the kinetic energy operator T in (8) we obtain 

(1 + /?) = ! [ a (a + /?+!) + 0+1] = ^, (14) 
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which in turn suggests that the von Roos ambiguity parameters are strictly 
determined (but not necessarily unique) as 



& = -\ , anda = 7 = -i. 
Hence, the (i-dimensional von Roos Hamiltonian reads 



H 



1 



m(q) 



dj + V(d) = - V d 



/ 1 



\m(q) 



(15) 



V d + V(d) (16) 



with 



v(a) = 



ldjdjm(q) 7 [djm(q) 



+ V(q). 



(17) 



4 m(q) 2 16 m(q) 3 

At this point, one may wish to subject such ambiguity parameters' settings 
(15) to Dutra's and Almeida's [11] reliability test on the exactly solvable one- 
dimensional Morse model (see equations (10)-(16) in [11]) . Such test shows that 
the ambiguous term v (a, /3,7,a) = ^/l/4 — 2q/c 2 = 1/4 (i.e., equation(16) in 
[11] for a = 7 = —1/4,(3 = —1/2) and classifies our ordering as a good-ordering 
(along with that of Zhu's and Kroemer's [18]. 



2.2 Corresponding d-dimensional radial Schrodinger Hamil- 
tonian 

We, in the forthcoming developments, shall assume the radial symmetrization 
of m (q) and V (q) in the <i-dimensional radially symmetric Schrodinger Hamil- 
tonian for d > 2. Under these settings, (17) and (18) imply 



where 



V(r) 



-d, 



m (r) 



dj + V(r) 



-V d 



m (r) 



■V d + V(r) (18) 



1 d 2 m(r) 7 [d r m(r)Y 



4 m(r) 2 16 m(r) s 



V(r); R9re(0,oo) 



(19) 



Recollecting that the (i-dimensional wave function for radially symmetric Schrodinger 
equation is given by 



* (r) - r- (d -^ 2 R rir4d (r) Y id , md {9 U 8 2 , 
would, in turn, when substituted in 

1 



m (r) 



V d + V(r) }0>(r) = E d t> (r) , 



(20) 



(21) 



results in the following <i-dimcnsional radial Schrodinger equation 



' d 2 £ d (£ d + 1) | m! (r) 
dr 2 r 2 m (r) 



2r 



d_ 

dr 



-m(r) V{r)-E d 



i d (r) = 0. 
(22) 
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Where id = ^+ {d — 3) /2 for d > 2, £ is the regular angular momentum quantum 
number, and ra r = 0, 1, 2, ■ ■ ■ is the radial quantum number. Of course, equation 
(23) is privileged with the inter-dimensional degeneracies associated with the 
isomorphism between angular momentum I and dimensionality d. On the other 
hand, moreover, the d = 1 (with l3r£ (0, oo) — > K 9 .t € (—00,00)) case 
can be obtained through the trivial substitutions id = — 1 and £<j = for even 
and odd parity, V = (— 1) respectively. Yet, a unique isomorphism exists 
between the 5-wave (£ = 0) energy spectrum in 3D and in ID. On this issue, 
the reader may wish to refer to, e.g., Mustafa and Znojil [32], and Mustafa and 
Mazharimousavi [13,14,31] and references cited therein. 



2.3 Corresponding PCT d-dimensional mapping 

In this section, we closely follow Mustafa's and Mazharimousavi's recipe dis- 
cussed in [14]. Where, a substitution of the form R{r) = m(r) 1 ^ 4 <p(Z (r)) in 
(22) would result in Z' (r) = \Jm{r), manifested by the requirement of a vanish- 
ing coefficient of the first-order derivative of <j> (Z (r)) ( hence a one-dimensional 
form of Schrodinger equation is achieved), and suggests the following point 
canonical transformation 

Z (r) = J T VMyjdy =f K T ,t d (Z (r)) = m (r)" 1/4 R n „j d (r) . (23) 



Which in effect implies 



d 2 | i d (i d + l) 
dZ 2 r 2 m (r) 



Veff (r) - E d } 4> nr ,e d (Z) = 0, 



(24) 



where 



V eff (r)=V(r)-Ud(r); U d (r) = ^ ( f 1} . (25) 

2rm (r) 

It should be noted, however, that the definition of Ud (r) in (25) is now more 
simplified than that in Eq. (8) of Mustafa and Mazharimousavi in [14]. 

On the other hand, an exactly solvable (including conditionally-exactly or 
quasi-exactly solvable) (i-dimensional time-independent radial Schrodinger wave 
equation (with a constant mass M (x) = m and h = 2m = 1 units) 

^ + Cd {C z d 2 + 1} +V(Z)- £ } ^ (Z) = (26) 

would form a reference for the exact solvability of the target equation (24) . That 
is, if the exact /conditionally-exact/quasi-exact solution (analytical/numerical) 
of (26) is known one can construct the exact /conditionally-exact/quasi-exact 
solution of (24) through the relation 

ij ^ + V(r)-Ud(r)-E d ^ ^^ + V(Z)-e, (27) 
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Where Cd is the (i-dimcnsional angular momentum quantum number of the 
reference Schrodinger equation. The reference -target map is therefore complete 
and an explicit correspondence (cf. e.g., Znojil and Levai [33] and Mustafa and 
Mazharimousavi [13,14]) between two bound state problems is obtained. 

A power-law position dependent mass of the form m (r) = sr v , for example, 
would imply a PCT function 

Z(r) = ^J r y v/2 dy=j^- ) r< v + 2 ^ =► t±V Z{r)=r ^r) (28 ) 
and hence equation (26) gives 



mr) = V ^= 2 ; v + -2 (29) 

2r 2 m(r) ( v + 2 ) 2 Z (r) 2 



Relation (28) in effect reads, 

A(A + 1) ( v , ^2 | ;: , A, (A,- I.) 

r 2 m (r) 

with 



Q + l) +V(r)-£a <=► +1J +V(Z)-e, (30) 



A = -i + |w + 2| 1 + l) 2 + A(, d (l d + 1) + 2v (1 - d); (31) 

However, for v = —2 =^> m (r) = ar -2 we obtain 

Z(r) = A^J t~ 1 dt = ^^r, (32) 



and hence 



U d {v = -2) = ed{ed + 1] -UAr,v = -2) 
£A£d + l) + d-l 



(33) 



Which would only add a constant to the left-hand-side of (28) to yield, with 
Cd = and/or Cd = — 1 (i.e., only s-states and/or d = 1 states are available 
from the right-hand-side of (28) ), 

V (r) + U d (v = -2) -E d ^ V (q) - e. (34) 



3 Concluding Remarks 

In this paper we have developed a Hermitian PDM-pscudo-momcntum opera- 
tor tlj = -i {F (q) dj + [djF (q)] /2} ,where F (q) = ±l/Jm (q). Hereby, the 
notion of PDM-pseudo-momentum operator is inspired by the fact that It has 
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an in-built regular momentum operator pj = —idj, which is recoverable at con- 
stant mass settings (i.e., M (q) = m =^> F (q) = ±1). Moreover, we have 
constructed our d-dimensional PDM-Hamiltonian, Hmm = n 2 + V (q) . 

On the other hand, upon intertwining our Hamiltonian, Hmm, with the von 
Roos d-dimensional PDM-Hamiltonian, H v r = Td(a,/3,j) + V (q), (cf.,e.g., 
Quesne [10]), we have observed that the so-called von Roos ambiguity parame- 
ters (i.e., a, (3 and 7) are strictly determined (i.e., (3 — —1/2 and a = 7 = —1/4), 
but not necessarily unique of course. Therefore, the von Roos d-dimensional 
PDM-Hamiltonian collapses into 

H vR H M m = -m (q)~ 1/4 d 3 m (q)~ 1/2 djiri (q)~ 1/4 + V (q) . (35) 

On the logistical supportive sides of our strict determination of the von Roos 
ambiguity parameters (3 = —1/2 and a = 7 = —1/4, we recollect that Bagchi 
et al [29], while analyzing the so-called quasi- free-particle problem, have used 
an intertwining relationship rjH = H\t] (where r\ is a Darbouxal first-order 
intertwining operator) and reported that such choices of the ambiguity parame- 
ters correspond to smooth mass functions m (a;) that signalled the formation of 
bound states. Yet, Koq et al. [30] have started with a = 7 = and (3 = — 1 with 
constant potential V (z) = V (equation (3) of Koc et al in [30]) in their study 
of transmission probabilities of the scattering problem through a square well 
potential with PDM barrier. However,they were forced to change the potential 
form (equation (4) of Koc et al in [30] ) into 

V{z) = V + -*L»{z)- 7 J£&) (36) 
8m {zf \ 4m (z) J 

which is exactly the same form of the effective potential that comes out from 
our eq.(17) with the new (3 = —1/2 and a — 7 = —1/4 parameters setting (of 
course one should mind the units used in this paper, h = 2m = 1). Moreover, 
Dutra's and Almeida's [11] reliability test resulted in classifying our ordering 
as a good-ordering (along with that of Zhu's and Kroemer's [18], and Li's and 
Kuhn's [19]). 

Therefore, the continuity conditions at the hetero junction boundaries and 
Dutra's and Almeida's [11] reliability test would ultimately single out Zhu and 
Kroemcr (a = 0, a = 7 = — 1/2, (3 = 0) [18] and our new ordering (/? = —1/2, 
a = 7 = — 1/4) as good orderings. 

On the least consequential research stimulant side, such ambiguity parame- 
ters' setting would, in effect, flourish a production-line for new sets of exactly- 
solvable, quasi-exactly solvable, and conditionally-exactly solvable target/new 
Hamiltonian models. The point canonical transformation (PCT) method used 
in this work exemplifies one of the methods that generate such spectrum of 
exact-solvability. For example, for a reference/old exactly-solvable 

V eff (Z) = ^±±A + v (Z) 
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in (30) there is a corresponding target/new exactly-solvable 



r 2 m (r) V2 / 



r 2 m (r) 

where A is given by (31) and v 7^ —2. Yet, a comprehensive number of illustra- 
tive examples on the generalized d-dimensional PCT is given by Mustafa and 
Mazharimousavi in [14]. Of course other methods designed to obtain exact- 
solvability do exist. Amongst, we may name the Lie algebraic method (cf., e.g., 
Levai in [29], intertwining operators related to supersymmetric quantum me- 
chanics (SUSYQM) method (cf., e.g., Quesnc in [10]), and the shape-invariance 
technique (cf., e.g., Quesne in [10] and Cooper et al. in [29]). 

On the feasible applicability side of our strictly determined von Roos ambi- 
guity parameters, the applicability of such ambiguity parameters' recipe should 
not only be restricted to Hcrmitian PDM Hamiltonians but also to a broader 
class of non-Hermitian PDM rj-weak-pseudo-Hermitian Hamiltonians (cf., e.g., 
Mustafa and Mazharimousavi [31] and related references cited therein). 
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